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ABSTRACT 

We estimate Cosmic Microwave Background (CMB) polarisation power spec- 
tra, and temperature-polarisation cross-spectra, from 7-ycar data of the Wilkinson 
Microwave Anisotropy Probe (WMAP). Foreground cleaning is implemented using 
needlet decompositions of sky maps for all WMAP channels, to produce maps for 
CMB temperature anisotropies (T-mode) and polarisation (E-mode and B-mode), for 
seven different years of observation. Power spectra are computed using a pixel- weighted 
scheme, from averages of independent cross-power estimates using foreground-cleaned 
maps for the different years. Error bars are estimated from the scatter of all possi- 
ble independent measurements of the band-averaged power spectra, and hence are 
independent of any detailed model of the WMAP noise. 

Our analysis technique yields a measurement of the EE spectrum on scales larger 
than about half a degree, that has smaller error bars than previous estimates and 
is in excellent agreement with theoretical expectations from the current concordance 
cosmological model. By comparison, the publicly available WMAP EE power spec- 
trum is higher on average (and significantly higher than the expected EE spectrum 
from the current concordance best fit cosmological model), although individual data 
points are compatible with ours within quoted errors. Our TE and TB measurements 
are also in good agreement overall with the WMAP ones, but our estimated error 
bars are smaller. Their compatibility with the theoretical expectations, however, is 
not quite as good on average as for EE, as some data points are off by a few stan- 
dard deviations. The upcoming Planck observations will help determine whether these 
data points are discrepant by reason of unmodelled residual systematics or foreground 
contamination, or indicate a tension in the best fit cosmological model. The EB and 
BB power spectra obtained in our analysis are compatible with zero, as expected for 
a standard cosmological model with low tensor to scalar ratio. 
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1 INTRODUCTION 

The Cosmic Microwave Background (CMB), relic radia- 
tion emitted when our Universe was about 380,000 years 
old, provides direct information about the origin and his- 
tory of cosmic structure. Hence, the current understand- 
ing of cosmological evolution is heavily based on observa- 
tions of the CMB. Over the last two decades, many ex- 
periments have accumulated observations on its tempera- 
ture anisotropies, imprints of the original density pertur- 
bations that later on gave rise to the large scale struc- 
tures observable today. Two very successful space missions, 
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the C osmic Background Explorer (COBE, iBennett et al.l 
1996), an d the Wilkinson M icrowave Anisotropy Probe 
(WMAP, IBennett et al1l2003h . complemented by several 
ground-based and balloon-borne experiments, such as CAT 
l|Baker etabll 19991). Boomeran g (jde Bernardis et alfcoool) 



Maxim a (lHananv et al. 200ol). A rcheop s dBenoit et al 
2003allbT). CBI (iPearson et al.ll2003l h VSA ilDickinson et al 



2004 



Rebolo et al. 



ACT l|Das et all 12011 



20041). ACBAR dReichardt et al 



a.b) 



and SPT (Kei sler et al 



2009), 



2011 



to mention just a few, have measured the CMB tempera- 
ture anisotropies at various angular scales and at various 
wavelengths. The CMB angular power spectra measured by 
these experiments agree remarkably well, and strongly sup- 
port a present 'concordance' cosmological model in which 
the Universe is spatially flat, with an energy density domi- 
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nated by the presence of about 70% dark energy and about 
30% matter, and in which structure forms from gravitational 
collapse of primordial adiabatic perturbations in the den- 
sity of the cosmological fluid. At present, this model does 
not require the additional admixture of primordial tensor 
perturbations (gravitational waves), although cosmological 
scenarios generically predict their existence. The observed 
CMB anisotropics are connected to the state of the cosmo- 
logical fluid at recombination in a way that depends on a 
small set of parameters specific of the cosmological model. 
CMB observations are of prime importance for constrain- 
ing these parameters of the model, as well as for testing its 
internal consistency. 

Temperature anisotropies alone, however, do not pro- 
vide the complete picture of the Universe. Independent 
information is needed to lift degeneracies between cos- 
mological parameter sets compatible with the CMB tem- 
perature power spectrum. Other cosmological probes such 
as the direct measure ment of the present expansion rate 
l|Freedman et al.l l200ll ), the observation of the expansion 
histo r y with supernova |Perlmutter et aLll 19991 ; lAstier et al.l 
120061 : iGuv et all |2010| ) , or baryonic oscillations t raced by 
the distribution of galaxies (Eise nstein et alj|2005l). provide 
comp lementary constraints ( Seliak. Slosar. fe McDonald! 
12009 ). 

On the side of the CMB itself, additional information 
are obtained by measuring CMB polarisation. The CMB 
is indeed partially polarised due to Thomson scattering of 
quadrupolar distributi on of photo ns on free electrons at the 
time of recombination (|Reesll968h , at a redshift z IBC of about 
1100. Large scale polarisation also arises from the scattering 
of CMB photons on free electrons after reionisation of the 
Universe at much lower redshift (z rc ion < 10). CMB polari- 
sation is an additional observable, which helps disentangling 
Sachs- Wolfe from Doppler contributions to the anisotropies, 
makes possible the estimation of individual contributions to 
the power spectra from scalar and tensor perturbations, and 
also helps constraining the epoch of reionisation. The ampli- 
tude of the polarisation, however, is significantly lower than 
that of temperature anisotropies, which makes its precise 
characterisation challenging. 

CMB polarisation was first detected at sub-degree 
angul ar scales by the DASI ground-based interferome- 
ter l|Leitch et all l2005h. It was subsequently measured 
by B oomerang jMontrqv et al.l |2006|: Piacentini et al.l 
|2006T I, Max ipol (TWu et al.ll2007F cBI (ISievers et alj|2007ft . 
CAPMAP jBischoff et al.ll200Sl). QUaD JPryke et alll2009l ; 



iBrown et al l 120091). BICEP dChiang et all I2010T ). QUIET 
lOUIET Collaboration et al.ll201ll ), and the Wilkinson Mi- 
crowave Anisotrop y Probe (WMAP, iNolta et all 120091 ; 
lLarson et afl|201lf ). 

The level of the CMB polarisation (at most a few per 
cent, depending on the angular scale), however, makes it 
easily contaminated by foreground emissions from both the 
Galaxy and extragalactic sources. Galactic emission, in par- 
ticular, can be polarised at the level of a few tens of per 
cent, i.e. ten times more than CMB anisotropies. The sig- 
nal to foreground ratio is thus less favourable for polarisa- 
tion than for intensity, in particular on large angular scales 
where Galactic synchrotron and dust emissions are strong. 
Contamination by foregrounds would result in CMB polar- 



isation power spectra larger than predicted by the current 
concordance cosmological model. 

Hence, contamination from polarised foreground signals 
has to be removed as much as possible for the accurate mea- 
surement of the temperature and polarisation angular power 
spectrum of the CMB. The effectiveness of a foreground 
cleaning technique is typically significantly improved by the 
use of prior knowledge of the emission properties of the con- 
taminant (such as the frequency dependence, typical angu- 
lar power spectrum of its emission, or availability of external 
templates). However, since the properties of foreground con- 
tamination are poorly known for polarisation (and, in partic- 
ular, no good template of polarised foregrounds is presently 
available), it is particularly crucial to develop and use data 
analysis tools that use only the minimum possible prior as- 
sumptions about foreground polarisation. 

Among the possible methods for CMB cleaning, the 
so-called Internal Linear Combination (ILC) method, first 
proposed for foreground cleaning i n the analysis of COBE- 
DMR data l|Bennett et al.l Il992h . and discussed subse- 
quently by many authors (see, e.g., iTegmarkl 1 19981 ; 
IDelabrouille fc Cardosol :2009), is a simple and effective 
way to combine multifrequency observations to extract 
the CMB while rejecting contamination by superimposed 
foreground signals. This method is based on two reason- 
ably safe assumptions. The first is that the amplitude of 
CMB emission is frequency independent in thermodynamic 
units, i.e. that the CMB emission law scales in frequency 
as the derivative of the blackbody spectrum of the cos- 
mic background. The second is that CMB fluctuations are 
not correlated to foreground signals. Under these assump- 
tions, the ILC estimates the CMB as a linear combina- 
tion of sky maps such that the variance of the estimate 
is minimum, while preserving unit response to the C MB 
(see, however, the appendix of IDelabrouille et al.l (|2009r i for 
a discussion of second order cor r ection s and biases, and 
iDick. Remazeilles. fc Delabrouille! (|2010l ) for a discussion of 
the impact of calibration errors). 

Component separation with an ILC method can be 
straightforwardly implemented either in real space or in 
harm oni c space jTegmark. de Oliveira-Costa. fc Hamilton 
20031 ; lEriksen et al l |2004 " Saha. Jain fc Soura dccp 



20061; 



Souradeep. Saha.. J air 
Saha. Prunet. Jain fc Souradeep j | 2008 



Souradeepj l201ll) Here, as in our previous work 



(IDelabrouille et all 120091 ; iBasak fc Delabrouille! 120121 1. 
we instead implement the ILC on a frame of spherical 
wavelets, called needlets. This special type of wavelets on 
the sphere provides good localisation in both pixel space 
and harmonic space because they have compact support in 
the harmonic domain , while still being very well localised 
in the pixel domain (|Narcowich. Petrushev fc Ward 2006; 
iMarinucci et alj 120081 ; iGuilloux. Fay, fc Cardosol I2009I '). 



Needlets have already been used in various analyses of 
WMAP data besides component separa tion and power 
spectr um estimation, for instance by IPietrobon et al 
(2008) to detect features in the CMB, and bv lRudiord et al 



(2009) to put limits on the non-gaussianity parameter /nl- 
Recently, we have used a needlet ILC on WMAP 7-year 
data to obtain a n estimate of the CMB temp erature angular 
power spectrum l|Basak fc Delabro uille 2012). In the present 
paper, as a natural extension to this work, we address the 
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problem of measuring as precisely as possible the CMB po- 
larisation power spectra from WMAP 7-year observations. 

The paper is organised as follows: In section [2] we de- 
scribe the methodology to estimate the CMB using an ILC 
on wavelet decompositions of multi- frequency sky maps. The 
implementation of this on WMAP 7-year data, and the re- 
sults for polarisation power spectra, are described in sec- 
tion [3] and section [4] We conclude in section [5] 



2 NEEDLET ILC ESTIMATE OF CMB 

A scalar field on the sphere such as CMB temperature 
anisotropies is conveniently expanded in usual spherical har- 



c can be modelled as, 



T(n) = 22 2~2 T lm Y lm (h) 



(1) 



The CMB polarisation field is usually specified using the 
Stokes parameters, Q and U, with respect to a particular 
choice of a coordinate system on the sky in relation to which 
the linear polarisation is defined. One can conveniently com- 
bine the Stokes parameters into the single complex quantity, 
P± = Q ± ill. Due to its rotation properties, one may ex- 

and P(h) in terms of spin-2 spherical harmonics, ±2Yi m (h) 

Goldberg et"ailll967h . as: 



p±m = Y, E p ^ 



lm ±2V; m (n 



(2) 



The description of CMB polarisation, however, tra- 
ditionally makes use of two scalar fields that are 
independent of how the coordinate system is ori- 
ented, and are related to the Stokes parameters by 
a non-local transformation (|Za ldarriaga fc Seliakl 1 19971 : 
iKamionkowski. Kosowskv. fc Stebbins! Il997l ). One of the 
fields, traditionally denoted as E, has even parity, whereas 
the other one, B, has odd parity (and hence is pseudo-scalar, 
rather than scalar). In harmonic space, the E and B modes 
of CMB polarisation are related to the complex polarisation 
fields P± as, 



P+2,lm = —{Elm + iBlm) 

and 



■(E h 



iBu 



(3) 



(4) 



Hence, one can fully characterise CMB anisotropies us- 
ing three Gaussianly distributed random scalar fields (T, E 
and B), without loss of information. In the following, maps 
of Q and U are systematically converted into maps of E 
and B by expansion onto spin-2 spherical harmonics, fol- 
lowed by an inverse spherical harmonic transform for Ei m 
and Bim independently. 



2.1 The CMB data model 

Denote X OBS ' c (h), (X = T,E,B) full-sky, multi-frequency 
temperature anisotropy and polarisation maps of the sky in 
n c different frequency bands (channels), such as those pro- 
vided by WMAP. The observed signals X OBS ' c (n) in channel 



X 



OBS,. 



(n)=/ dn A , b c (h.n) X slG ' c (h') + X N - c (h) (5) 



where X slG ' c (h) is the signal (sky) component, itself decom- 
posed in the sum of CMB and foreground components, 



X (n) = a X (n) + X (n), 



(6) 



a c being the CMB calibration coefficient for the channel c. 
Up to calibration uncertainties, a c — 1 for all WMAP chan- 
nels. If, in addition to WMAP data, we use external data 
sets which serve as foreground templates to help foreground 
subtraction, as done in the present work, the coefficients a c 
vanish for such data sets (i.e. the ancillary maps contain no 
CMB anisotropies). 

The beam function b c (h.fi), represents the smoothing 
of the signal due to the finite resolution of the observations. 
Assuming for simplicity that the beams are circularly sym- 
metric (a good approximation for WMAP data), b c (n.h') 
depends only on the angle 9 — cos^ 1 (n.n') between the di- 
rections n and n', and can be expanded in terms of Legendre 
polynomials, 



b c ( 



°° 21 + 1 



(T) 



The last term, X N,c (n), in equation \E\ represents the 
detector noise in channel c, and is not affected by the beam 
function. For a spherically symmetric beam, equation can 
be recast straightforwardly in the spherical harmonic repre- 
sentation, as: 



X 



OBS,. 



c , c V CMB . , c v FG,c . v N,c 

-a Mta + b i A i m + Ji i m 



(8) 



where Xi m stands for the three modes Ti m , Ei m and Bi m of 
temperature and polarisation in harmonic space. 



2.2 Implementation of the needlet transform 

Considering that each channel observes the sky at a differ- 
ent resolution, the maps are first convolved /deconvolved, in 
harmonic space, to the same resolution: 



Xlm 



h 
b? 



X; 



OBS., 



(9) 



Each of these maps Xf m is then decomposed into a set of 
filtered maps X^ represented by the spherical harmonic 
coefficients, 



h-i X [m , 



(10) 



where the filters h\ , serving for localisation in the harmonic 
space, are chosen in such a way that 



E 



(ii) 



The reconstruction of the original maps Xf m from the collec- 
tion of the filtered maps Xf£, representing each a different 
scale, is performed using the same set of filters. In terms of 
h\, the spherical needlets are defined as, 



*ifc(n) = \A^E E hlY lm (h)Yi m (£ jk ), 



(12) 



1=0 m=-l 
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where denote a set of cubature points on the sphere for 

scale j. In practice, we identify these points with the pixel 
centr es in the HEALPix pixelisation scheme (jGorski et al.l 
2005). Each index k corresponds to a particular HEALPix 
pixel, at a resolution parameter nside(j) specific to that 
scale j. The cubature weights Ajk are inversely proportional 
to the number Nj of pixels used for the needlet decomposi- 
tion, i.e. Xjk = 4r- The needlet coefficients for CMB fields 
X(h) are denoted as, 



Pfk = 



X(h) *jfe(n) dVlf, 



Xjk ^E ^E hj bi Xi m Y; m (£ 



jk) 



(13) 



1=0 m=-l 



The linearity of the needlet decomposition implies that the 
needlet coefficients /3j k corresponding to the filtered map 
obtained from the harmonic coefficients are a linear 

combination of the needlet coefficients of individual compo- 
nents and noise at HEALPix grid points £jk'- 

a X,c c oCMB . oFG.c . a N ,c 

Pjk = a Pjk + Pjk + Pjk 

where, 



(14) 



oCMB 
Pjk 



jk 



Pjk 



^ fe E E hibiX™ B Yi m (£j k ) 

1=0 m=—l 

'max I 



1=0 m=-l 
imax I 



A ^E E H$xZ?Yim(&k) 

1 = m=-l 1 



(15) 



2.3 Implementation of the needlet ILC 

The ILC estimate of needlet coefficients of the cleaned map 
is obtained as a linearly weighted sum of the needlet coeffi- 



cients 



'j k ! 



and of 



where indices c and c' of elements of the matrix ^jfe 1 ] 
the vectors Qjk and a are written down explicitly for clarity. 
More compactly, we have: 



UJjk = 



[%'] 



[%■] 



(19) 



where Qjk is the vector of ILC weights to be applied to the 
needled coefficients of all input observations at scale j and 
in pixel k, a is the CMB 'mixing vector' (a vector of n c 
entries all equal to unity for inputs in thermodynamic tem- 
perature), and ^Jj; 1 ] an estimate of the inverse covariance 
of the needlet coefficients of the n c observations at pixel k 
of scale j. 

The NILC estimate of the cleaned CMB needlet coeffi- 
cients is: 



P. 



NILC 



= Pjk + 2^ "Jk [Pjk + Pjk 



(20) 



The elements of the covariance matrix for scale j at 
pixel k, R^ k = (^PjkPjk^i are obtained each as an average 
of the product of the relevant computed needlet coefficients 
over some space domain T>k centred at k. In practice, they 
are computed as 



Rx,jk = — E w * ( fc > k ') Pfk C Pfk 



X,c 



(21) 



where the weights Wj(k,k') define the domain X\. A sensi- 
ble choice is for instance Wj(k, k') — 1 for k' closer to k than 
some limit angle, and Wj(k,k') = elsewhere, or alterna- 
tively, Wj(k,k') shaped as a Gaussian beam of some given 
size that depends on the scale j (which is what we do here). 

Finally, the NILC estimate of the cleaned CMB map 
can be reconstructed from cleaned CMB needlet coefficients 
using the same set of filters that was used to decompose the 
original maps into their needlet coefficients. The NILC CMB 
temperature map is then 



X NILC (h) 



(22) 



oNILC c a 

Pjk = 2^ U J k P. 



jk 



(16) 



where ojj k is the needlet weight for scale j and frequency 
channel c, at the pixel k of the HEALPix representation of 
the needlet coefficients for that scale. Under the assump- 
tion of decorrelation between CMB and foregrounds, and 
between CMB and noise, the empirical variance of the error 
is minimum when the empirical variance of the ILC map 
itself is minimum. The condition for preserving the CMB 
signal during the cleaning is encoded as the constraint: 



(17) 



The resulting needlet ILC weights Qj k that minimise the 
variance of the reconstructed CMB, subject to the constraint 
that the CMB is preserved, are expressed as: 



u jk 



(18) 



with 

V NIL' 



, V CMB . yRFG . yRN 



(23) 



where the harmonic coefficients residual foreground (Xf^ ) 
and residual noise (Xf^) are given by: 

AL FG = E E Pf™ h i Y m tefc) (24) 



and 



X™ = E E \f^jk Pfk h\ Ylmfek). 



(25) 



Equations [221 and [231 imply that the NILC estimate of CMB 
contains some residual foreground and noise contamination. 



3 WMAP 7- YEAR NEEDLET ILC MAP 

The WMAP satellite has observed the sky in five frequency 
bands denoted K, Ka, Q, V and W, centred at 23, 33, 41, 
61 and 94 GHz respectively. After 7 years of observation, 
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Table 1. List of needlet bands used in the present analysis. 



Band index 








nside 


1 








50 


32 


2 





50 


100 


64 


3 


50 


100 


150 


128 


4 


100 


150 


250 


128 


5 


150 


250 


350 


256 


6 


250 


350 


550 


512 


7 


350 


550 


650 


512 


8 


550 


650 


800 


512 


9 


650 


800 


1000 


512 



the released data includes temperature anisotropy and po- 
larisation maps obtained with ten difference assemblies, for 
7 individual years. One map is available, per year, for each 
of the K and Ka bands, two for the Q band, two for the V 
band and four for the W band. These sky maps are sampled 
using the HEALPix pixelisation scheme at a resolution level 
(nside = 512), corresponding to approximately 3 million sky 
pixels. 

We work on band-averaged maps of T, E and B for 
the five frequency bands, complemented, for temperature 
only, by thre e foreground templates (dust at 100 m icrons, as 
obtained by ISchlegel. Finkbeiner. fc Davis! dl998h . the 408 
MHz synchrotron map of Haslam et alj (1 1981 ), an d the com- 
posite all-sky H-alpha map of lFinkbei ncr (2003)). All sky 
maps are convolved /deconvolved in harmonic space, to a 
common beam resolution (corresponding to that of the W 
frequency channel of WMAP 7- year release). 

Each of these maps is then decomposed into a set of 
needlet coefficients. For each scale j, needlet coefficients of 
a given map are stored in the format of a single HEALPix 
map at degraded resolution. The filters h\ used to compute 
filtered maps are shaped as follows: 



for P < / < P 

iul l min *» 1 ^ 'peak' 



for / = I 



peak , 



for P <T I < P 



For each scale j, the filter has compact support between the 
multipoles P min and P max with a peak at P peak (see figure [T] 
and table[T]). The needlet coefficients f3f k are computed from 
these filtered maps on HEALPix grid points £j k with resolu- 
tion parameter nside equal to the smallest power of 2 larger 
than P max /2. The estimates of needlet coefficients covari- 
ance matrices, for each scale j, are computed by smoothing 
maps of products of needlet coefficient /3| fc /3| fe with Gaus- 
sian beams. In this way, an estimate of needlet covariances 
at each point k is obtained as a local, weighted average of 
nearby needlet coefficient products. The full width at half 
maximum (FWHM) of each of the Gaussian windows used 
for this purpose is chosen to ensure the computation of the 
statistics by averaging about 1200 samples or more, result- 
ing from a trade-off between the localisation of the estimates 
(which requires small windows) , and the accuracy of the es- 




200 400 600 800 1000 

Multipole(l) 

Figure 1. Needlet bands used in the present analysis. The black 
line shows the normalisation of the needlet bands, i.e. the total 
filter applied to the original map after needlet decomposition and 
synthesis of the output map from needlet coefficients. 



timate (which require large windows). Choosing a smaller 
FWHM results in inaccuracy in the covariance estimates, 
and hence ILC bias. Choosing a larger FWHM results in 
less localisation, and hence loss of effectiveness of the needlet 
approach. 

Using these covariance matrices, ILC weights are com- 
puted for each of T, E and B, for each scale j and for each 
pixel k of the needlet representation at scale j. For each of 
T, E and B, a full sky CMB map, at the resolution of the 
WMAP W channel, is synthesised from the NILC needlet 
coefficients. 



4 WMAP 7- YEAR NEEDLET ILC SPECTRUM 

As already mentioned above, the foreground-cleaned maps 
of T, E and B obtained in this way are not fully exempt 
from contamination by residual foregrounds and noise. The 
recovered CMB map (T, Q and U) at 60 arcmin resolution 
is displayed in figure [21 Residual foreground contamination, 
albeit small, is visible along a narrow strip on the Galac- 
tic plane of the maps. Noise contamination is seen from the 
larger variance of the polarisation maps away from the eclip- 
tic poles. 



4.1 Minimisation of the impact of noise 

Noise biasing in our estimated power spectra is avoided by 
producing, for each of T, E and B, an independent CMB 
map for each of the 7 individual years of observation, and 
computing the CMB power spectra exclusively from cross- 
products between maps from different years. Each data point 
in our spectra is thus obtained as an average of 21 possible 
cross-year spectra. All maps, however, are obtained using 
the same set of needlet weights, determined using the co- 
added 7 year observations. 



4.2 Minimisation of the impact of foregrounds 

Unlike the residual noise, the residual foreground emission 
in each of 7 cleaned maps is the same, as each map observes 
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Figure 3. Apodisation kernel for Galactic mask used in our anal- 
ysis 



■ 0.030 mKCMB 



We use Mioo 10 o, Mis°,ib°, M2o°,2o° to compute auto and 
cross angular power spectra involving T, E and B field of 
CMB respectively. 

Finally, for power spectrum estimation, small gaps in- 
duced by the point source mask are filled-in by an interpola- 
tion procedure using the values of CMB anisotropics in the 
neighbouring unmasked pixels. 




0.030 mKCMB 



Figure 2. The NILC estimate of temperature anisotropies(top) 
and, Stokes parameters Q (middle) and U (bottom) of CMB ob- 
tained by implementing the NILC on the 7-year band average 
maps, at nside = 512 and degraded to 60 arcmin resolution. 



the same sky emission, and is produced using the same linear 
combination for all years. 

For safety, we use a conservative mask for estimat- 
ing CMB angular power spectra, that combines the point 
source mask provided by the WMAP collaboration with an 
apodised Galactic mask, applied directly on the CMB maps 
after the needlet ILC. We choose this option (rather than 
masking before the ILC) so that we produce full sky CMB 
maps that can also be used for other purposes than power 
spectrum estimation. 

The Galactic mask Mg g used in practice is: 



Mg r ,g w {SA) = ' 







for 0° ^ 



I for 



for > Or + 1 



4.3 Impact of calibration errors 

An important assumption of the ILC is that the frequency 
scaling of the CMB is known. However, calibration coeffi- 
cients for each channel, which are a multiplicative factor for 
each frequency, introduce an uncertainty in the frequency 
scalings of t he CMB component in the presence of c alibra- 
tion errors (|Dick, Remazeilles. fc Delabrouillel koioh . This 
effect is particularly strong in the high signal to noise ra- 
tio regime. Considering the relatively low signal to noise of 
WMAP polarised maps, this issue can safely be ignored here. 

Beam uncertainties induce similar biases as calibration 
uncertainties, except that these biases are scale dependent. 
Here again, such biases are not the main source of error in 
our final polarisation spectra, as their impact is small in 
comparison to the uncertainties due to instrumental noise 
for polarisation measurements with WMAP. 

These issues connected to the exact response of the de- 
tectors, however, will require specific attention with upcom- 
ing more sensitive observations of the polarised CMB, if our 
method is to be used for the analysis of these future data 
sets. 



4.4 Noise-weighting 

Residual noise in our CMB maps is inhomogeneous, primar- 
ily because of non-uniform sky coverage, with higher number 
of observations in the directions of ecliptic poles and rings 
at 45° ecliptic latitude. 

The estimation of band-averages of the CMB power 
spectrum (in bins of I), results in lowering the effective reso- 
lution in harmonic space (as compared to computing power 
spectra for every single value of I). This results in the capac- 
ity to measure the bin-averaged power spectrum with some 
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Table 2. Comparison of our estimate of binned angular power 
spectrum of _E-mode of CMB polarisation with that provided by 
WMAP team. The quantities tabulated are D t = l(l + \)Cf B /2-K 
and AD, = 1(1 + l)ACf ;E /2n. 
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localisation in pixel space, and hence compute the band av- 
eraged power as weighted averages of local estimates. We 
take advantage of this opportunity and follow the general 
idea of the noise-weighting over needlet domains already 
performed in our previous an alysis of WMAP intensity maps 
IjBasak fc Delabrouillell2012j). also investigat ed for tempera- 
ture in a method paper bv lFav et al.l i|200Sft . 

The contribution of instrumental noise, in each bin of I, 
at each location of the sky, is estimated using 100 indepen- 
dent realisations of the reconstruction of the CMB by the 
NILC. This local noise variance is then used to compute a 
weighted average of the band-averaged power, estimated lo- 
cally from the covariance of maps filtered with the spectral 
window of interest. This method yields more precise esti- 
mates of CMB power spectra in those angular scales where 
residual noise dominates over residual foreground signal, a 
feature of particular interest for polarisation measurement 
on small scales, for which the noise is the main source of 
error. 

Note that the apodised Galactic masks used in the 
power spectrum estimation (see equation ([26} and figure [3}, 
also act as effective weights over the the pixel space (with 
vanishing weights at the lowest Galactic latitudes). 

We are using here the me thodology already described 
in lBasak fc Delabrouillel l|2012h . The adaptation for polari- 
sation and cross spectra is straightforward, and we refer to 
our previous work for implementation details. 



4.5 Results 

We now present the estimated polarisation power spectra 
we obtain on the basis of WMAP 7 year observations. The 
error bars in our estimates include the total statistical error 
of the estimator (noise and cosmic variance terms) . 
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Figure 4. The red circles show the angular power spectrum of 
E-mode of CMB polarisation estimated using 7 years of observa- 
tions of WMAP. The dark blue squares show the 7-year angular 
power spectrum of _B-mode of CMB polarisation published by the 
WMAP collaboration. The solid black line shows the theoretical 
angular power spectrum for WMAP best-fit A-CDM model. The 
top panel uses a linear scale in the horizontal axis, and the bottom 
panel a logarithmic scale. 



4-5.1 EE angular power spectrum 

Figure [4] shows the estimated auto-angular power spectrum 
for the _E-mode of CMB polarisation. Our result for i?-mode 
CMB power is slightly higher than that obtained by WMAP 
collaboratiorQ in the first bin (on very large angular scales, 
corresponding to I = 4), although the two are compatible 
within error bars (see table • Considering that we use a 
different mask, different methodology, and different weight- 
ing scheme, the small difference between the two measure- 
ments is not too surprising. In most of the rest of the scales, 
our measurements are lower than the WMAP result, and 
in better agreement with the theoretical expectations (as- 
suming the concordance cosmological model is correct). The 
difference between our measurement and that of the WMAP 
team may point out that the latter could be biased by excess 
residual foreground contamination. 



1 http://lambda.gsfc.nasa.gov/data/map/dr4/dcp/spectra/wmap_ee_spectrum_7yr 
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Figure 5. The red circles show the cross-angular power spectrum 
of CMB temperature anisotropics (T-mode) and £-mode of CMB 
polarisation estimated using 7 years of observation of WMAP. 
The dark blue squares show the 7-year angular power spectrum 
of U-mode of CMB polarisation published by the WMAP col- 
laboration. The solid line shows the theoretical angular power 
spectrum for WMAP best-fit 'concordance' A-CDM model. The 
top panel uses a linear scale in the horizontal axis, and the bottom 
panel a logarithmic scale. 



4-5.2 TE cross-spectrum 

Figure [5] shows the estimated cross-angular power spectrum 
for temperature anisotropy of CMB and i?-mode of CMB 
polarisation. The power spectrum obtained using our anal- 
ysis agrees well with that provided by the WMAP collabo- 
ration and is in general consistent with the ACDM model 
for WMAP best-fit cosmological parameters. Our error bars, 
computed from the internal scatter of 21 individual measure- 
ments in each bin, are smaller than those provided by the 
WMAP collaboration^. We note, however, that four of our 
data points between I — 50 and I = 120 (at I = 56, 65, 98, 
111) are significantly above the WMAP best fit model (up 
to ~ 4 standard deviations away) (see table[3}. Whether this 
indicates that the WMAP best fit model is not a good fit 
to the observations, or that our error bars are impacted by 
un-modelled systematics, is at present not clear, ft will be 
investigated further with upcoming data sets, such as those 
from the Planck mission. 



Table 3. Comparison of our estimate of binned cross-angular 
power spectrum of CMB temperature anisotropics (T-mode) and 
E-mode of CMB polarisation with that provided by WMAP 
team. The quantities tabulated are D' t = lCf E /2n and ADj 
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We also note that our cross-spectrum estimate seems 
to be consistently slightly lower on average than that of the 
WMAP collaboration, as well as the theoretical expectation 
from the concordance model, around the peak at I — 300. 
Individual points agree within error bars, but there seems 
to be a systematic difference that remains to be explained. 



4-5.3 TB cross-spectrum 

Our result for the TB cross-spectrum, compared to the 
WMAP collaboration one0 is shown in figure [5] Theoreti- 
cally, this cross-spectrum is supposed to vanish (to preserve 
the parity symmetry), and significant departure from zero 
would be the sign of either unknown systematics in the mea- 
surement (including residual foregrounds), or new physics. 

Although our estimated TB cross-power spectrum is in 
general in good agreement with the expectation of a vanish- 
ing cross-power, a few of the data points shows significant 



http://lambda.gsfc. nasa.gov/data/map/dr4/dcp/spectra/wmap_ tc_spb iittrjn(^irabq^ 
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Figure 6. The red circles show the cross-angular power spectrum 
of CMB temperature anisotropics (T-modc) and B-mode of CMB 
polarisation estimated using 7 years of observation of WMAP. 
The dark blue squares show the 7-year angular power spectrum 
of _E-mode of CMB polarisation published by the WMAP collab- 
oration. The top panel uses a linear scale in the horizontal axis, 
and the bottom panel a logarithmic scale. 



E 

I 

rxr-~ 



E 



o 



-5x10 



-1x10 



1x10 



5x10 



-5x10 



-1x10 



■ WMAP7yr 

• NILC7yr 



5x10 - 



100 200 300 400 500 600 700 
Multipole(l) 



■ WMAP7yr 

• NILC7yr 



10 100 
Multipole(l) 

Figure 7. The red circles show the cross-angular power spectrum 
of _B-mode and B-mode of CMB polarisation estimated using 7 
years of observations of WMAP. The top panel uses a linear scale 
in the horizontal axis, and the bottom panel a logarithmic scale. 
We use the same legend and color code for data points in this 
figures as in Figs. \E\ [6] and [8] although WMAP data points are 
absent here. 



departure from zero (few sigma) (see tabled}. This discrep- 
ancy between model and measurement, although small, is 
presently not fully understood, although a possible expla- 
nation is effects of beam im perfections, inducing a leak- 
age of T into B (see, e.g., iKaplan fc Delabrouille! l2002i ; 
iHu. Hedman. fc Zaldarriagal 120051 ; iRosset et all 120071 . for 
discussions of such effects) . A detailed investigation of these 
discrepancies requires very accurate modelling of the mea- 
surement and of the data processing pipeline, and is beyond 
the scope of the present paper. 



4-5.4 EB cross-spectrum 

The WMAP collaboration has not provided cross-angular 
power spectra for E and B-modes of CMB polarisation. Our 
estimated EB cross-power spectrum, displayed in figure [7] 
and tabulated in tableO is compatible with zero as expected 
from the current cosmological best-fit model. 



4-5.5 BB angular power spectrum 

Finally, the power spectrum of CMB B-modes is of special 
interest, as it provides one of the most promising means of 
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Figure 8. The red circles show the angular power spectrum of 
B-mode of CMB polarisation estimated using 7 years of observa- 
tions of WMAP. The dark blue squares show the 7-year angular 
power spectrum of i?-mode of CMB polarisation published by the 
WMAP collaboration. 
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Table 4. Comparison of our estimate of binned cross-angular 
power spectrum of CMB temperature anisotropics (T-mode) and 
B-mode of CMB polarisation with that provided by WMAP 
team. The quantities tabulated are = lCf B /2-K and AflJ 



Table 5. Our estimate of binned cross-angular power spectrum 
of _B-mode and _B-mode of CMB polarisation. The quantities tab- 
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detecting primordial tensor modes in the early universe, and 
constrain models of inflation. Significant effort is currently 
undertaken to prepare the measurement of B-modes with 
a future space mission. One such mission, the Cosmic Ori- 
gins Explorer (COrE), has been proposed to ESA within 
Cosm ic Vision 2015-2025 (|The COrE Collaboration et all 
l201ll ). Missions with similar objectives, but different designs, 
have been proposed to NASA (see, e.g., iBaumann et al.l 
2009). Contamination by foregrounds is one of the main 
worries for this measurement, and the investigation of the 
severeness of the contamination, as well as the development 
and validation of component separation methods adapted to 
the challenge of measuring CMB _B-modes, has focused sig- 
nifica nt attention recently (iTucci et al.l 120051; IStivoli etal 



2006: 



lAmblard. Coorav. fc Kaplinghatl 120071; jBetoule et al 
lEfstathiou. Gratton. &i Paci 



2009; Dunkl ev et all 20091 
20091 ; IStivoli et al.ll2010h . 

Clearly, the WMAP mission lacks the sensitivity to 
place a meaningful limit on the B-mode polarisation. As 
expected, neither our result nor that of the WMAP collab- 



Table 6. Comparison of our estimate of binned angular power 
spectrum of B-mode of CMB polarisation, with that provided by 
WMAP team. The quantities tabulated are D t = 1(1 + l)C^ B /2n 
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5.700e-08 


3(5 


-9.478e-08 


1.437e-07 


1.055e-07 


1.248e-07 


74 


-4.331e-08 


7.960e-08 


2.272e-07 


2.560e-07 


124 


-4.248e-07 


2.177c-07 


4.432e-07 


5.966e-07 


174 


-4.622e-07 


-1.519e-06 


8.142e-07 


1.165e-06 


224 


1.858e-07 


4.507e-07 


1.594e-06 


2.092e-06 


274 


-3.411e-06 


-2.052e-06 


2.811e-06 


3.581e-06 


324 


-2.517e-06 


-1.631e-06 


4.979e-06 


5.952e-06 


374 


9.446e-06 


7.744e-06 


9.846e-06 


9.643e-06 
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oratiorQ do show any detection of B-modes (see figure [8] 
and table [6]). The method used in the present paper, how- 
ever, has been also tested on realistic simulations of future 
observations with COrE. On such simulations, it allows to 
reject foreground contamination effectively enough to mea- 
sure tensor to scalar ratio of 10 -3 , limited by the sensitivity 
of the observations rather than foreground contamination. 



5 CONCLUSIONS: 

Following our previous work on WMAP temperature maps, 
we have computed CMB power spectra for polarised WMAP 
observations. CMB polarisation maps for T, E, and B are 
obtained from WMAP observations using linear combina- 
tions that minimize the variance of the recovered CMB on 
wavelet (needlet) domains (the Needlet ILC method), that 
are subsequently used to compute CMB power spectra. 

Our analysis differs substantially from that of the 
WMAP team: we use all WMAP channels, use a needlet 
ILC over the full range of harmonic modes, and produce 7 
independent maps for each of T, E and B, from the different 
years of observation. Error bars do not rely on a model of the 
WMAP noise, but instead are computed directly from the 
internal scatter of all possible (21 in our case) independent 
band-averaged cross-spectra. 

We find that our EE power spectrum is in excellent 
agreement with the expectations from the current 'concor- 
dance' cosmological model, while the 7-year WMAP spec- 
trum publicly available on the Lambda web sit^fl seems to 
be systematically higher. Our statistical error bars on EE 
are lower comparable to those of the published WMAP re- 
sult in most of the multipole bins under consideration (and 
are, actually, 10-20% lower on average, except for the first 
bin). On TE and TB, our error bars are also smaller than 
the WMAP ones, but we note that a few data points deviate 
by few sigma from the expected best-fit value, a feature that 
will need to be investigated with future data sets. Finally, 
our measurements of the EB and BB spectra are compati- 
ble with zero. 
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